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Предложен способ оценивания линейного фазового сдвига ОРйМ сигнала, основанный на методе наименьших квадратов и его 
модификации. Проведена оценка эффективности данного способа. Рассмотрена возможность его реализации на микросхемах 
программируемой логики. 


Введение 

Формирование ОРЭМ (ОііИо<юпа1 Ргециепсу 
Оіѵіиоп Миіііріехкщ - ортогональное частотное 
мультиплексирование) сигнала происходит путем 
обратного преобразования Фурье от поднесущих, 
модулированных по закону квадратурной ампли- 
тудной модуляции. Ріа выходе обратного преобра- 
зования Фурье-сигнал периодический. 

Промежуток времени, для которого амплитуда и 
фаза поднесущих остаётся постоянной, называют 
ОРОМ символом. Для правильного декодирования 
сигнала длительность символа должна быть больше 
или равна периоду сигнала на выходе обратного пре- 
образования Фурье. Для компенсации межсимволь- 
ных искажений, вызванных многолучевой природой 
распространения сигнала, в начало символа поме- 
щается циклический префикс ( сусііс ргефіх [1]), пред- 
ставляющий собой копию конца символа, рис. 1. 



На рис. 1 Г РРТ - период сигнала полученного пу- 
тем обратного образования Фурье (РРТ - РаМ Рои- 
гіег ТгашГогт), Г ср - длительность циклического 
префикса, Г 5утЬо1 - длительность ОРОМ символа. 

На приемной стороне символ оцифровывается 
на интервале определения сигнала, равного дли- 
тельности 7рр Т . В результате прямого дискретного 
преобразования Фурье для каждой поднесущей 
ОРОМ символа определяется амплитуда и фаза. 


Если дискретная выборка принятого символа не 
совпадает по времени с Т ш , то на выходе декодера 
каждая поднесущая ОРОМ символа получит фазо- 
вый сдвиг, что приведет к ошибке декодирования. 

Компенсация фазового сдвига производится 
СОКОІС (СОогсІіпаіе Коіаііоп ОІфіаІ СопгрЩег) 
процессором [2]. Для корректной работы СОКОІС 
процессора необходимо провести оценку фазового 
сдвига каждой из поднесущих или временного 
сдвига ОРОМ символа. 

В данной статье предложен способ оценки вре- 
менного сдвига ОРОМ символа и, как следствие, 
оценки линейного фазового сдвига поднесущих 
ОРОМ сигнала. 

Оценка временного сдвига 

При временном сдвиге ОРОМ символа фаза 
поднесущих изменяется следующим образом: 

С к =РРТ{8„), 

: 2 л кт 

С к е М = РРТ(8 п _ т ), 

где РРТ— оператор прямого быстрого преобразова- 
ния Фурье, - дискретная комплексная выборка 
из ОРЕ>М символа, С к - комплексные точки, соот- 
ветствующие созвездиям квадратурной амплитуд- 
ной модуляции (поднесущие), N — количество то- 
чек преобразования, количество поднесущих, т - 
смещение начала выборки, к - номер поднесущей, 
п — номер выборочного отсчета ОРОМ символа. 

Для оценивания т необходимо использовать 
пилотные поднесущие. На этих поднесущих пере- 
дается определенная последовательность, которая 
известна на приемной стороне. Приемная сторона 
может сравнивать параметры принятых пилотных 
поднесущих с эталоном. 


148 



Технические науки 



Рис. 2. Плотности распределения фазы суммы гармонического сигнала и шума: — по (3), - - - по ( 4) 


Для нахождения оптимального алгоритма оце- 
нивания т необходимо знать плотность распреде- 
ления зашумленного сигнала. 

Предположим, что на каждую пилотную подне- 
сущую воздействует аддитивный узкополосный 
нормальный стационарный шум 

1(0 = Д0со8[й) 0 Г-Д(0], (1) 

где А(і) - огибающая шума, щ - эквивалентная ча- 
стота шума, 3(1) - фаза шума, і— время. 

Каждую поднесущую можно выразить как гар- 
монический сигнал 

•КО = К С 08 (® с О, 

где А т - амплитуда сигнала, со с - частота сигнала. 

Тогда сумму сигнала и шума можно представить 
в виде: 

5(0 + КО = А т со8(®Д) + Д0со8[<зу -Д(0] = 

= V (0 со8(ш 0 ? — у/(0) = Ѵ(і)соѣ(ю с 1-Ѳ(і)), 

где Аа>=со-(о 0 , Ѳ(1)=Аа>1+у/(і), Ѵ(1) - огибающая 
суммы сигнала и шума, 9(1) — фаза суммы сигнала 
и шума. 

Согласно [3] совместная плотность вероятности 
для огибающей и фазы суммы гармонического сиг- 
нала и шума (1) в один и тот же момент времени 
равна: 

/(Г,ѳ) = 

= -^ехр[--і у (К 2 + АІ- 2 ѴА т со8(0))], 

2ко 2(7 

ѵ>0, (2) 

где ст 2 - дисперсия шума ( 1 ). 

Для нормального случайного процесса с симме- 
тричной спектральной плотностью огибающая А(1) 
и фаза (р(і+т) независимы как в совпадающие, так 
и в разные моменты времени [3]. Исходя из этого, 
при оценке т необходимо использовать только фа- 
зу пилотных поднесущих. 


Проинтегрировав выражение (2) по всем воз- 
можным значениям Кот 0 до со, получим плотность 
распределения фазы суммы гармонического сигна- 
ла и шума [3] 

№ = 

1 “Ій 2 I 7С а 2 С05 2 (Ѳ) 

= -^— е 2 [1+ ѵ2жасо8(Ѳ)Ф(о-со8(Ѳ))е- ], 


-к <Ѳ <л, 


( 3 ) 


где а - отношение сигнал/шум, а = — Ф(ф) - ин- 

<7 

теграл вероятности, который определяется как: 


Ф( 2 ) 




Ах. 


Выражение (3) можно упростить для малых и 
больших значений а. 

Так, если й=0, то плотность распределения фа- 
зы становится равномерной. 

Когда а>>1, с некоторым приближением мож- 
но считать 

ехр(-ідг) » 0, со8(0) «1, 8Іп( 0) я Ѳ, Ф(а) «1, 


тогда выражение (3) может быть преобразовано: 


т 



\ а 2 Ѳ 2 

2 


а » 1. 


( 4 ) 


Выражение (4) соответствует нормальной плот- 
ности распределения. На рис. 2 изображены плот- 
ности распределения фазы суммы гармонического 
сигнала и шума при разном отношении сигнал/шум. 

Зная плотность распределения фазы суммы гар- 
монического сигнала и шума, можно выбрать ме- 
тод оценивания параметра т. В качестве оценки 
параметра т допустимо использовать оценку мак- 
симального правдоподобия [4]. 

Предположим, что шум для разных пилотных 
поднесущих некоррелирован. Тогда оценивание 
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параметра т можем быть проведено на основе ме- 
тода наименьших квадратов. 

Подход с помощью метода наименьших квадра- 
тов к задаче оценивания содержится в теореме 
Гаусса. Она утверждает, что если ошибки I, некор- 
релированы и имеют нулевое среднее значение 
Е[І ]= 0 и одинаковую дисперсию Е[2 2 ]=с г 2 , то оп- 
тимальными выборочными оценками параметров 
являются значения, минимизирующие сумму ква- 
дратов расхождений между наблюдаемыми значе- 
ниями и подбираемой моделью 


Н(8 1 ,8 2 ,-:,Е к ) X, ( У, 8 х Х п 8 2 Х і2 ... 8 к Х ік ) , (5) 

і=1 

где 8і - подбираемые параметры, минимизирую- 
щие Н(е ь 8 2 ,—,8 к ), Ь - количество отсчётов функ- 
ции, по которым производится оценка. 

Для рассматриваемого случая все условия со- 
блюдаются. Тогда (5) можно переписать: 

Н(т) = X (Ф**. - той(1 80)) 2 , (6) 

і=і 



Рис. 4. Вероятность неправильного оценивания пп в зависимости от отношения сигнал/шум при 256-точечном преобразова- 
нии Фурье 
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где N - количество точек преобразования Фурье, 
і - номера пилотных поднесущих, I - количество 
пилотных поднесущих, <р ікст - фазы принятых пи- 
лотных поднесущих. 

В выражении (6) вычисляемое значение фазы 
должно находится в пределах ±180°, так как на при- 
ёмной стороне происходит измерение главного (в 
пределах ±180°) значения фазы. 

Для оценивания параметра т необходимо по- 
добрать такое т, чтобы Н(т) была минимальной. 
Стандартный подход к решению этой задачи за- 
ключается в дифференцировании Н{т) по т и на- 
хождении такого т, при котором получившееся 
выражение равно нулю. Подробно метод и резуль- 
таты описаны в [4] . 

Для нашего случая такой метод не подходит, так как 
имеется множество локальных минимумов (рис. 3). 

Находить параметр т можно перебором всех 
возможных его значений. 

При реализации вычисления выражения (6) 
при помощи микросхем программируемой логики 
желательно избежать применения операции возве- 
дения в квадрат, т. к. такие процедуры требователь- 
ны к аппаратному ресурсу. Эту операцию можно 
заменить операцией взятия модуля. Поэтому вме- 
сто выражения (6) можно использовать 

т> я) = X ( Р.^е - (=Ц М ) ПНХК180) (7) 

і=1 ' /Ѵ 

При воздействии шума на сигнал абсолютный 
минимум (рис. 3) становится менее выраженным, а 
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в некоторых случаях и вообще исчезает. В связи с 
этим при определении т могут возникать ошибки. 
На рис. 4 представлены графики вероятности при- 
нятия ошибочного решения при вычислении оце- 
нок по выражениям (6) и (7). 

Заключение 

Оценивание линейного фазового сдвига ОРОМ 
сигнала можно производить при помощи метода 
наименьших квадратов. При этом необходимо на- 
ходить минимум суммы квадратов отклонений фа- 
зы пилотных поднесущих от идеальных значений. 
Более простым, с точки зрения вычислений, явля- 
ется нахождение минимума суммы модулей откло- 
нений фазы пилотных поднесущих от идеальных 
значений. Если данный способ оценивания реали- 
зован на микросхеме программируемой логики, то 
существенно экономится аппаратный ресурс ми- 
кросхемы. Так, операция возведения в квадрат 8- 
разрядного числа занимает 30 логических элемен- 
тов, а операция взятия модуля такого же числа за- 
нимает 8 логических элементов. Для 16-разрядных 
чисел возведение в квадрат и взятие модуля зани- 
мают соответственно 126 и 16 логических элемен- 
тов. То есть замена операции возведения в квадрат 
на операцию взятия модуля позволяет уменьшить 
объем вычислений для 8 -разрядных чисел более 
чем в 3 раза, а для 16-разрядных - более чем в 7 раз. 

При замене операции возведения в квадрат на 
операцию взятия модуля для одной и той же веро- 
ятности принятия ошибочного решения необходи- 
мо увеличение отношения сигнал/шум на 0,15 дБ. 
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